Dzyaloshinskii— Moriya interaction: How to measure its sign in weak ferromagnetics? 
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Three experimental techniques sensitive to the sign of the Dzyaloshinskii-Moriya interaction are 
discussed: neutron diffraction, Mossbauer 7-ray diffraction, and resonant x-ray scattering. Classical 
examples of hematite (a-Fe2 0a) and MnC03 crystals are considered in detail. 
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Weak ferromagnetism (WF) of antiferromagnetics is a 
classical example of an initially small and controversial 
physical problem that later produces a strong impact on 
the general picture of magnetic phenomena. From the 
very beginning, the modern theoretical consideration of 
WF developed by Dzyaloshinskii and Moriya was based 
on symmetry arguments, both phenomenological [l|, Q 
and microscopic 0, Q. It was shown that appropriate 
crystal symmetry allows the following term in the 
interaction of two antiferromagnetic sublattices Si and 
S 2 



D-[SixS 2 ], 



(1) 



which favors to (usually small) canting angle between 
Si and S2', here D is a vector parameter of the 
Dzyaloshinskii-Moriya interaction. Possible directions 
of D were found [J] for different local symmetries. 
Significant progress was recently achieved in ab initio 
calculations of D (see [5| and references therein). 
Besides fundamental interest, the Dzyaloshinskii- 
Moriya interaction is a very important ingredient of 
magnetoelectric effects with possible applications to 
spintronics. 

The canted spin arrangement is just responsible for 
WF. Both the magnitude of WF and the canting angle 
are proportional to \D\ and therefore it seems that the 
sign of D is not important at all. According to Eq. 
(UJ), the sign obviously depends on our choice which of 
sublattices is 1 or 2 and therefore it is usually claimed 
that the sign is conventional. From the phenomenological 
point of view this is true because in macroscopic 
theory the phase of antiferromagnetic arrangement is 
not fixed relative to the crystal lattice. However at 
the atomic level the phase can be fixed owing to the 
Dzyaloshinskii-Moriya interaction and the sign of this 
interaction is crucial for relation between the local crystal 
structure and magnetic ordering. For instance, this sign 
determines the handedness of spin helix in crystals with 
the noncentrosymmetrical B20 structure [il 0] ■ In this 
paper we show how one can measure it in classical WF 
crystals like cv-Fe2 03 or M11CO3. 



Let us rewrite Eq. JT|) in a more invariant form not 
depending on any arbitrary choice of sublattices. If two 
atoms with spins si and s 2 are located at the points n 
and r 2 , then we can add the following scalar to the energy 
of their interaction 



Tjk m sijS 2 k(ri - r 2 ) 



(2) 



where an antisymmetric tensor, Tjkm = ~Tkj m , 
characterizes interaction of spins Si and s 2 through 
intermediate crystal matter. The properties (in particular 
symmetry) of the intermediate matter determines the 
properties of this tensor including symmetry restrictions 
on its tensor components. It is well known also that 
any third-rank antisymmetric tensor is equivalent to 
a second-rank pseudo-tensor A nm : Tjkm — tjk n A nm 
where tjkn a unitary antisymmetric pseudo-tensor (A nm 
changes its sign under inversion). The relation between 
D, A nm and T jkm is given by 



i(ri - r 2 ) r 



(3) 
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Using the well known symmetry restrictions on the 
third-rank antisymmetric tensors [8| we can obtain from 
Eq. ([2]) all the symmetry restrictions on D found in 
In particular, A nm = 0, Tjkm = and D m — if the 
points ri and r 2 are related by inversion symmetry (rule 
1 from [H). If there is an n-fold rotation axis (n > 2) 
along Ti — r 2 then D is parallel to T\ — r 2 (rule 5 from 
1), etc 

However, there is an important principal difference 
between D and Tjkm'- tensor Tjk m (or pseudo-tensor 
A nm ) can be considered as a field on the lattice, it should 
be invariant relative to all the symmetry operations of the 
space group. In particular, it is determined by the same 
parameters at all equivalent lattice points; of course, one 
should take into account corresponding crystallographic 
operations connecting those equivalent points: rotations 
(changing orientations of the principal axes) and space 
inversions (changing signs of all components of Tjkm)- 
On the contrary, the pattern of vector D on the 
lattice cannot be obtained by pure crystallographic 
operations and some additional consideration is needed 
(see discussion of La2Cu04 in \9L 

We conclude this short introduction with a remark 
that Eq. @ cannot be used for quantitative description 



of WF; modern first-principles theoretical considerations 
are more appropriate [5|. Nevertheless this expression can 
be used for better understanding of symmetry aspects of 
the problem and now we will show that this is really the 
case. 

The appearance of the antisymmetric third-rank tensor 
suggests an idea that there is some chiral effect behind 
this. In fact, it is known for a long time [Il[ that 
the Dzyaloshinskii-Moriya interaction can produce long- 
period magnetic spiral structures in ferromagnetic and 
antiferromagnetic crystals lacking inversion symmetry. 
This effect was suggested for MnSi and other crystals 
with B20 structure [a, 0| and it has been carefully proved 
that the sign of the Dzyaloshinskii-Moriya interaction, 
hence the sign of the spin helix, is determined by 
the crystal handedness |12j. More delicate situation 
with chirality occurs in typical WF crystals with 
R3c symmetry (a-Fe203, M11CO3, etc.) which are 
centrosymmetric. 

At first let us consider classical WFs, carbonates 
of transition metals, for instance MnCC>3 [l3j. In its 
primitive rhombohedral unit cell with the space group 
i?3c, there are two Mn atoms at crystallographically 
equivalent inversion centers 3, vertices (0, 0, 0) and body- 
centers (i, i, i). Atoms at the vertices and body-centers 
have almost opposite magnetic moments lying in the 
planes normal to the threefold axis. According to the 
crystal symmetry the moments should not be exactly 
opposite and in fact the moments are slightly canted 
so that the resulting WF moment is along one of three 
twofold axes. There are also two carbon atoms at points 
and (|, |, |). 

The physical origin of WF in M11CO3 is a weak 
relativistic interaction between spins in the lattice with 
the i?3c space group. But what is the structural origin 
of this i?3c symmetry? If we consider only Mn atoms, 
the symmetry of lattice would be R3m and WF would 
be impossible. This symmetry will not change if the 
carbon atoms are taken into account. Only oxygen atoms 
change i?3m symmetry to i?3c; thus their configuration 
is crucial for the value (and sign) of the Dzyaloshinskii- 
Moriya interaction and it is worthy of a more careful 
consideration. 

In MnCC>3, there are six hexagonal Mn layers per the 
lattice period along the threefold axis, so that atoms of 
the next layer is just under the centers of triangles formed 
by atoms from the previous layer, and layer sequence is 
ABCABC..., like in the fee lattice. These equidistant 
layers has z-coordinates equal to z = 0,g,^,|,|,l; in 
each layer all the spins are parallel and are lying in the 
layer plane. The spin of neighboring layers are almost 
opposite. Here and below the standard hexagonal setting 
of the rhombohedral lattice is used |14j . 

Considering the first two layers one can see that 
between them, at z = , there is a layer of oxygen and 
carbon atoms; the point symmetry of this layer is 32, 
just because of low symmetry of oxygen positions: oxygen 
atoms are at the 18e positions with point symmetry 
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right-hand twist left-hand twist 

Phc. 1: Right-hand and left-hand twists of moments between 
layers alternating along z axis. Triangles indicate the threefold 
axis normal to the figure plane; all possible directions of 
WF moments in MnCC>3 (twofold axes) are shown by small 
arrows. Big arrows are spin directions in neighboring layers 
at different z levels for the case when external magnetic field 
is applied in horizontal direction and DM is positive. A bold 
small arrow indicates the direction of WF moment. 

2 and with coordinates equivalent to (xo,0, 4) where 
xo ~ 0.27. Carbon atoms are at the 6a positions (0, 0, 4) 
with point symmetry 32. 

It is very important that this intermediate layer is 
noncentrosymmetric and therefore both vector D and 
tensor Tjk m can have some nonzero values. For pairwise 
interaction of spins from different layers, tensor Tjk m 
has symmetry I, but being averaged over all pairs it 
has of course the symmetry of the intermediate layer. 
For symmetry 32, tensor Tjk m is determined by two 
independent parameters, say T yzx and T xyz but only 
the latter leads, according to Eq. @, to a twist angle 
between the Mn spins lying in the first and second layers. 
The sign of this twist angle is just the sign of T xyz . This 
twist violates the right-left symmetry and, in a figurative 
sense, we can say that the intermediate C-0 layer is 
chiral. 

The next intermediate C-0 layer is at z ~ j, i.e. 
between the Mn layers at z = i and z — i, and for this 
layer all the components of tensor Tjkm change sign due 
to inversion centers at z — i. Thus the "chirality" of this 
layer is opposite to that of the first intermediate layer, 
hence the small twist angle between spins is also opposite 
and for the Mn layer at z = ^ the spin orientation 
exactly coincides with the spin orientation for the Mn 
layer at z — (see Fig. []} . And then this repeats from 
layer to layer. We see that two alternating local twists, 
left and right, between neighboring Mn layers result in 
macroscopic canting angle between magnetic sublattices 
in centrosymmetric WF crystals. 

If we change the sign of xo so that x' = 1 — xo ~ 
0.73, then the layer "chirality" changes to opposite, T xyz 
change the sign and the twist angle also change the sign. 
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From pure crystallographic point of view both values, 
xo and x' Q , are equivalent, they simply correspond to 
the lattice origin shifted at a half-period, from (0,0,0) to 
(0,0,^). In the crystallographic databases, both values 
of xo are cited for a-Fe203; nevertheless sometimes the 
first-principle results for WF in a-Fe2C>3 do not indicate 
which value of xo they really adopt (see for instance 
[H[). We see also that in R3c crystals the sign of the 
Dzyaloshinskii-Moriya interaction changes to opposite at 
one half-period, therefore the idea to measure this sign 
by the Mossbauer absorption cannot be correct. 

According to i?3c symmetry there are six possible 
orientations of the WF moments (along plus and 
minus directions of three twofold axes; in Fi g. [H 
they are shown by short arrows, like in |14|). 
Application of an external magnetic field along one of 
these directions makes the corresponding ferromagnetic 
domain energetically favorable. And if the orientation 
of the ferromagnetic domain is fixed then the 
Dzyaloshinskii-Moriya interaction fixes the phase of 
antiferromagnetic sequence of moments in this domain 

Now we are ready to consider the main item of this 
paper: how to measure the sign of the Dzyaloshinskii- 
Moriya interaction in WFs? First of all, a strong 
enough magnetic field should be applied to obtain 
the single domain state where the Dzyaloshinskii- 
Moriya interaction pins antiferromagnetic ordering to the 
crystal lattice. Next, single-crystal diffraction methods 
sensitive both to oxygen coordinates and to the phase 
of antiferromagnetic ordering should be used. In other 
words, one should observe those Bragg reflections hk£ 
where interference between magnetic scattering on Mn 
atoms and nonmagnetic scattering on oxygen atoms is 
significant. There are three suitable techniques: neutron 
diffraction, Mossbauer 7-ray diffraction, and resonant x- 
ray scattering. We will discuss now their advantages and 
disadvantages. 

Because of the layered magnetic structure alternating 
along z-axis, the reflections with strong magnetic 
scattering correspond to the reciprocal lattice vectors 
Hhu with odd I. At first let us consider scattering on 
oxygen atoms; the expressions for the oxygen structure 
amplitudes F ox (H) are similar for all three techniques 
and looks like (i = 2n + 1): 

F ox (H) = 2A ox [cos 2ir(hx ox + 1/4) + cos 2ir(kx ox +1/4) 

+ cos 2iv{hx ox + kx ox — 1/4) 
— 8A ox (—l) n+1 sin 7r/i2; oa ; sin nkx ox smn{h + k)x ox ,{4) 

where A ox is proportional to the oxygen atomic 
scattering factor (for x-rays) or to the nuclear scattering 
length (for neutrons); it is practically real because oxygen 
absorption is very small for thermal neutrons or hard 
x-rays. There is no contribution from carbon fl4j . Here 
and below some factors (such as atomic and magnetic 
formfactors, the Debye-Waller factor, etc) are omitted 
because corresponding expressions are well known and 



implemented into routine computer programs used for 
diffraction experiments. It is clear from this equation that 
one should measure reflections with hk(h + k) ^ 0. 

In the case of neutron diffraction, one can adopt 
the standard technique using the polarization ratio 
R, i.e. the ratio of reflection intensities for incoming 
neutrons with spin <x parallel and antiparallel to the 
direction of applied magnetic field. For I — 2n + 1 
this ratio is given by the following expression containing 
interference between nuclear scattering by oxygen atoms 
and magnetic scattering by Mn atoms 



R(hkt) = 



\F ox {H) + a ■ Q{H) 
\F ox (H) - a ■ Q(H) 



\F ox (H) - s DM A mag k{2h + k)\ 2 
\F ox {H) + s DM A mag k(2h + k)\< 



(5) 
(G) 



where Q(H) is the magnetic structure amplitude for 
reflection H, Q(H) cx M H - H(H ■ M H )/H 2 , M H 
is the correspondent Fourier harmonic of the vector 
field describing the electron-magnetization distribution, 
sum is the sign of the Dzyaloshinskii-Moriya interaction 
between the first two layers of Mn atoms and A mag 
includes all routine factors describing neutron magnetic 
scattering. Calculating (JB]) from ([5]) we took into 
account the geometry shown in Fig. [T] (i.e. <x is 
directed horizontally, Mh is directed vertically so that 
er ■ Mh = if small canting is neglected, etc.). 
In particular, factor sdm appears just because the 
phase of antiferromagnetic sequence is fixed by the 
Dzyaloshinskii-Moriya interaction. All other factors in 
Eq. ^ are more or less known and one can determine 
sdm from rather rough measurements giving R(hk£) < 
1 or R(hk£) > 1. Notice that there is an additional 
condition for this measurements, 2h + k =/= 0, which 
appears in Eq. ^ from H ■ Mh ^ 0. To avoid 
confusion we should notice that expression (QJ obeys 
threefold symmetry whereas magnetic scattering does 
not, because of the external field applied perpendicular 
to the threefold axis. 

The technique of polarized neutron diffraction was 
used for measurements of the s ign of small angular 
deviations of moments in MnF2 [18| but in that case 
the deviation is introduced by the single-spin anisotropy 
[HI]. To the best of our knowledge, there were no 
attempts to measure with this technique the sign of the 
Dzyaloshinskii-Moriya interaction. 

In hematite (a-Fe203), the situation is slightly more 
complicated because iron atoms are in the positions 
(0, 0, zp e ) with the point symmetry 3 and neighboring 
Fe layers are coupled either ferromagnetically or 
antiferromagnetically. Like in MnC03, the i?3c symmetry 
is induced by oxygen atoms. The Fe layers coupled 
feromagnetically are related by inversion and, according 
to the Morya rules, there is no canting between them 
(Tjk m = for inversion centers). The antiferomagnetic 
neighboring layers interact via an oxygen layer with 
symmetry 32 and alternating right and left twists of their 
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moments lead to a macroscopic WF moment. Eq. © 
transforms to 



R{hk£) 



SDAiA mag k(2h + k) cos 2ir£z 



Fe 



+ SDMA mag k(2h + k) cos2irizFe\ 2 



(7) 



The additional factor cos27r£z^ e (where zp e k, 0.355) 
allows us to change the value and sign of magnetic 
scattering just changing £. A suitable reflection (210 
rhombohedral, i.e. 213 hexagonal) had been studied in 
[l7^ ; however, as it was noted in that paper, the result 
was controversial: in the great majority of observations, 
i?(213) = l/i?(213) whereas i?(213) = i?(213) would 
be expected from the symmetry of hematite. Thus we 
cant extract the sought sign of the Dzyaloshinskii-Moriya 
interaction and more careful experiments are needed. 

It should be again emphasized that our symmetry- 
based arguments are only qualitative: the ab initio 
calculations for a-Fe2C>3 show that the torque induced 
by neighboring antiferomagnetic layer is opposite to the 
total torque [5j. This confirms importance of suggested 
direct measurements of the sign of the total torque. 

The Mossbauer diffraction can be used in a similar 
way. In this case, there is no need to vary the photon 
polarization and one can study intensity Itki{E) of 
reflections as a function of photon energy E: 



Ihke(E) = \F ox (H) — SDMB rnag (H ,E)\ 



(8) 



where B mag (H , E) is the magnetic Mossbauer structure 
factor of the hk£ reflections with £ = 2n+l, hk(h+k) ^ 0. 



The function B mag (H , E) is well known |20j and its 
real part, which interferes with the first term in ([5]), 
changes sign when E passes through resonances provided 
by the hyperfine splitting of nuclear levels. This should 
facilitate the observation of interference between the 
two terms. For 57 Fe both terms in may be of the 
same order of magnitude. The Mossbauer diffraction was 
observed in many crystals including WF a-Fe2C>3 (2l| 
and FeBC>3 [22j (an analog of M11CO3) but these studies 
were concentrated mainly on pure magnetic scattering 
rather then on its interference with scattering on oxygen 
atoms. Contrary to neutrons, the Mossbauer diffraction 
can be used for very thin layers but the number of 
possible crystals is rather limited by the list of suitable 
Mossbauer isotopes. 

Another promising approaches to the sign 
measurements can be related with resonant x-ray 
diffraction, i.e. diffraction near x-ray absorption edges. 
It is sensitive both to structural and magnetic ordering 
especially near L absorption edges (see recent surveys 
[23H25j |y However, in d magnetic metals, the K edge 
is the only appropriate for diffraction, and for this 
edge magnetic scattering is several orders of magnitude 
smaller than conventional charge scattering by electrons. 
Thus, for reflections of I = 2n + 1, hk(h + k) ^ type, 



the intensity is given by eq. ((5]) with the second term 
much smaller then the first one. Therefore the reliable 
observation of interference between two terms will be 
perhaps very difficult. 

However, resonant x-ray scattering provides another 
nontrivial way to measure sdm ■ The asymmetric oxygen 
environment of transition metals induces some additional 
anisotropy of their non-magnetic scattering amplitude, 
so that just owing to this anisotropy the reflections with 
£ = 2n + 1 can be excited even if hk(h + k) = 0. These 
reflections do not exist out of the resonant region and 
they are referred to as "forbidden reflections" . There is no 
direct contribution to forbidden reflections from oxygen 
atoms, equation Q gives zero, but the sign of the induced 
anisotropy depends on asymmetrical arrangement of the 
oxygen atoms and, correspondingly, the sign of the non- 
magnetic structure amplitude of forbidden reflections 
with £ — 2n + l,hk(h + k) = is proportional to 
the sign of xo- For crystals with i?3c symmetry these 
reflections were first observed in a-Fe203 [2(|. Then it 
was predicted that there should be some "chiral" dipole- 
quadrupole contribution to these reflections (2?} and 
interference between different contributions (including 
magnetic scattering) have been studied in detail for 
a-Fe203 and Cr2 03 crystals [28j|. It was shown that 
the azimuthal dependence of reflection intensity could 
be strongly influenced by this interference (especially 
for the weak 009 reflection, see Figs. 11 and 12 from 
[1^1 ) and orientation of antiferromagnetic moment was 
determined in a-Fe2 03 from the observed azimuthal 
dependence (without external magnetic field). Exactly 
the same measurements in orienting magnetic field would 
allow us to determine the sign of the Dzyaloshinskii- 
Moriya interaction. 

The only problem with the last method is that we 
should rely on the sign of the x-ray anisotropy of iron 
atoms calculated with rather sophisticated computer 
codes. However it was proved experimentally for Ge (29| 
that those codes (we used FDMNES Q) are rather 
reliable. It is worth noting that resonant x-ray diffraction 
and Mossbauer diffraction are element sensitive and 
moreover the former can distinguish orbital and spin 
contributions to magnetic moments. 

In conclusion, we see that the experiments similar 
to those needed for the sign measurements had been 
already performed (some of them long time ago) for all 
three considered techniques. Therefore we believe that 
this paper will stimulate these measurements in different 
types of the weak ferromagnetics. 
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